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We calculate the far-from-equilibrium dynamics and thermalization both for the quantum 
and the classical 0(iV)-model. The early and late-time behavior can be described from 
the 2P/-loop expansion for weak couplings or the nonperturbative 2PI-1/N expansion 
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of the effective action beyond leading order. A comparison with exact simulations in 1 + 1 
dimensions in the classical limit shows that the 2PI-1/N expansion at next-to-leading 



order gives quantitatively precise results already for moderate values of N. We derive a 



criterion for the validity of the classical approximation and verify it by comparing far-from- 



equilibrium quantum and classical dynamics. At late times one observes the expected 
deviations due to the difference between classical and quantum thermal equilibrium. 

r£3 ' In recent years we have witnessed an enormous increase of interest in the dynamics of 
quantum fields out of equilibrium. Strong motivation comes from a wide range of applica- 
tions including current and upcoming relativistic heavy-ion collision experiments, phase 
transitions in the early universe or the dynamics of Bose-Einstein condensation. Directly 
simulating quantum fields in real time, such as solving the Schrodinger equation for the 
wave functional is prohibitively difficult and one has to find suitable approximations. The 



nonequilibrium time evolution is inherently nonperturbative in the sense that approxima- 



tions based on a finite order in standard perturbation theory break down at sufficiently 
late times. Practicable approximations for nonequilibrium dynamics may be based on the 
two-particle irreducible (2PI) generating functional for Green's functions Recently, 
the 2PI effective action has been solved for a 1+1 dimensional scalar quantum field theory 
at next-to- leading order in the 2P/-loop expansion |3j and in the 2PI-1/N expansion [|J. 
Both the far-from-equilibrium early-time behavior and the late-time physics of thermal- 
ization were successfully described. For a recent review on the use of the 2PI effective 
action in nonequilibrium field theory see Ref. f|. 

A unique possibility to calculate the exact time evolution, which includes all orders in 
loops or 1/N, is provided by the classical statistical field theory limit. The exact evo- 
lution (up to statistical errors) of correlation functions can be constructed by numerical 
integration and sampling of initial conditions from a given probability distribution func- 
tion. On the level of correlation functions classical and quantum evolution equations are 
remarkably similar^ and the same approximation schemes and initial conditions can be 
applied. This aspect of classical field theory has been stressed in Refs. [pHnj. In Ref. || 
this is applied to the next-to-leading order classical 4 -model. In Ref. [|nj it is shown that 



1 Of course, in contrast to the quantum theory the classical limit suffers from Rayleigh- Jeans divergences 
and has to be regulated. In 1 + 1 dimensions such divergences are absent in 00-correlation functions ij. 
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the 1/N expansion at next-to-leading order converges to the exact result by increasing 
N already for moderate values of N. Apart from benchmarking approximation schemes 
employed in quantum field theory, the classical field limit is of great practical importance 
and often applied for the approximate description of nonequilibrium quantum fields. In 
this note I want to elaborate on Ref. ||10|| , done together with G. Aarts, and study the 
conditions under which far-from-equilibrium quantum dynamics can be reliably described 
by classical fields. 

We consider a real iV-component scalar quantum field theory with a \(<f) a (j) a ) 2 /(4\N) 
interaction in the symmetric phase (a = 1, . . . ,N). There are two linearly independent 
two-point functions which can be related to the anti-commutator and commutator of two 
field operators [p~T|.[T| 

F ab (x,y) = ([(f) a (x),(j) b (y)} + }/2 , p a6 (x, y) = i([(/>a(x), Mv)]-) i 1 ) 

Here F is the "symmetric" propagator and p denotes the spectral function. The classical 
equivalent of the spectral function is obtained by replacing the commutator by the Poisson 
bracket. For the analytic presentation we consider here the three- loop expansion of the 
2PI effective action for N = 1, which has been employed in Ref. || to study late-time 
thermalization in a quantum field theory. Numerical results from the 1/N expansion of 
the IP I effective action at next-to-leading order for N > 1 are shown below. For spatially 
homogeneous fields the dynamics of the Fourier transformed F and p is described by [|TT|,f4f 

dt z E p (t x ,t z ;p)F(t z ,t y ;p) 



[dl + M 2 (Q] F(t x ,t y ]p) 



+ 



[dl+M 2 (t x )} p(t x ,t y ;p) 



dt z E F (t x ,t z ; p)p(t z ,t y ] p) 

■x 

dt z E p (t x ,t z ;p)p(t z ,t y ;p) 



(2) 



which are exact for known E_p, E p . From the three- loop 2PI effective action the effective 



mass term is M 2 (t x ) = m 

\ 2 I" dq dk 



Sf^I, ty] P) 
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+ (A/2) J F(t x , t x ; q) and the self energies are 
F(t x ,t y ;p- q-k) 



{27i) d (27r) d 



F(t X , ty] q)F(t X , ty] It) - - p{t X , ty] q)p(t x , ty] k) 



^p(txi ty, p) 
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dq dk 

(27i) d (2vr) d 



p(t x , t y ;p-q-k) 



(3) 



F(t X , ty] q)F(t X , ty] k) - — P(t x , ty] C^) p{t X , ty] k) 



The classical statistical field theory limit of a scalar quantum field theory has been studied 
extensively in the literature. An analysis along the lines of Refs. |1^JT3| J9|JI0|] shows that 
all equations (|2|)-(|3|) remain the same in the classical limit except for differing expressions 
for the self energy: 

A 2 f dq dk 



^F{tx, ty] p) 



classical limit 



6 J (27r) d (27r) 



F(t x , t y ;p-q- k)F(t x , t y ] q)F(t x , t y ] k), (4) 
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Figure 1. Nonequilibrium classical time evolution from the three-loop approximation 
of the 2PI effective action for A/m 2 = 1. Shown is the mode temperature T'(t;q), 
as defined in the text, with T'(t 



0;q)/m = 5 | 



The equivalent approximation 
for the corresponding quantum field theory has been used in Ref. fU to demonstrate 
the late-time approach to quantum thermal equilibrium. In contrast, the classical field 
theory approaches classical thermal equilibrium as one clearly observes from the approach 
T'{t; q) — > T c \ which corresponds to classical equipartition. 



^pij^Xl ty, p) 



classical limit 
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<iq dk 

(2vr) d (2vr) 



1 p(t x ,t y ;p- q-k)F(t x ,t y ;q)F(t x ,t y ;k). (5) 



One observes that the classical self energies are obtained from the expressions in the 
quantum theory by dropping terms with two spectral (p) components compared to two 
statistical (F) functions. In particular, it becomes obvious that the leading order equa- 
tions (similarly for leading-order large N, Hartree or mean field) are identical for the 
quantum and the classical theory, and the inclusion of direct scattering effects is crucial. 

The solution of the above classical evolution equations is shown in Fig. [I] for Gaussian 
initial correlations with temperature T'(t = 0;q)/m = 5 and A/m 2 = 1.0 (For the 
numerical implementation see Ref. M.) The "mode temperature" is defined by |7| 



T>(t xlP )=d tx d ty F(t x ,t y ,p)\ t 



(6) 



One observes that the system relaxes to a final temperature T c \/m = 5.5. The thermal- 
ization time turns out to be very large and we find an exponential late-time relaxation to 
thermal equilibrium with rate ^( therm ) ~2x 10~ 4 for T(t; q = 0) ||14j| . We emphasize that 
the exponential behavior with similarly long thermalization times are found as well from 
exact simulations as in Ref. J7|Q. 



2 For the classical theory we employ a lattice regularization with spatial lattice spacing ma s = 0.4 corre- 
sponding to a fixed momentum cutoff A = n/a s . We observe that at sufficiently late times the contribu- 
tions from early times to the dynamics are effectively suppressed. This fact is shown in detail in Ref. p4| 
and has been employed in Fig. |l| to reach the very late times. 
3 Note that these authors employ a fixed lattice cutoff n/a s — Airm. 
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Figure 2. Left: far-from-equilibrium evolution of the two-point function F(t,t;p) for 
various momenta p in units of tur. The coupling is X/6N = 0.5 m R for N = 10. One 
observes a good agreement between the exact (dashed) and the next-to-leading order 
classical result (full) ||10|| . The quantum evolution is shown with dotted lines for momenta 
p < 2p ts , for which the classicality condition (|j) is approximately valid. Right: A very 
sensitive quantity to study deviations is the time dependent inverse slope T(t,p) defined in 
the text. When quantum thermal equilibrium with a Bose-Einstein distributed particle 
number is approached all modes get equal T(t,p) = T qm , as can be observed to high 
accuracy for the quantum evolution PUlOfl. For classical thermal equilibrium the defined 
slope remains momentum dependent and T'(t,p) becomes constant (cf. Eq. @). 



The classical field approximation is expected to become a reliable description for the 
quantum theory if the number of field quanta in each mode is sufficiently high. Comparing 
Eqs. (^) and (Q) one observes the sufficient condition for 

3 

classical evolution: F(t x , t y ; q)F(t x , t y ; k) > - p{t x , t y ; q)p(t Xl t y ; k) . (7) 

To obtain an estimate on the occupation numbers we define a time-dependent effective 
particle number n(t;p) and mode energy e p (t) [11]] 

n(t x ;p) + -= (F(t x ,t y -p)d tx d ty F(t x ,t y ]p)) , e p (t x ) = — ^— r — 

I V /\t x =t v \ !<{t x ,t y ;p) J \ tx=ty 

With these definitions F(t x , t x ; p) = (n(t x ; p) + l/2)/(e p (t x )). Note that p = for t x = t y 
and at unequal times we find \p(t x , t y ; p)| < max[(e p (t))~ t ^ tx t J. The latter relation 
has indeed to be valid for the free theory or if the weakly coupled quasiparticle picture 
applies, which is underlying the above particle number and mode energy definitions. Time 
averaged over an oscillation period 2 r n/e p {t) one can translate (0) into an estimate on the 
lower bound for the effective particle number: 



2 



3 

> - or n(t;p) > 0.37 . (8) 
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This limit agrees rather well with what is found in thermal equilibrium. For a Bose- 
Einstein distributed particle number = (e €//T — with temperature T one finds 
ny(e = T) = 0.58, below which deviations from the classical thermal distribution become 
sizeable. 

A similar estimate can also be obtained beyond the weak coupling regime from the 
1/N expansion of the 2PI effective action at next-to-leading order. In Fig. |2] we consider 
the full next-to- leading order time evolution in the (2PI-) 1/N expansion @,[nj. The 
far-from-equilibrium initial particle number, n (p), is described by a Gaussian n (p) = 
*4exp(— 7^2-(H — Pts) 2 ) peaked around p ts = 2.5?71r with o = 0.25m^ and A = 4. We add 
to no a thermal "background" ut with initial temperature T = 4m#. Here tur is the 
one-loop renormalized mass in vacuum (n = 0). 

For the current initial distribution the classicality condition (^j) is approximately valid 
for momenta p < 2p ts with n(t = 0;p = 2p ts ) ~ 0.35 (n(t = 0;p = p ts ) ~ 4.5) and a 
slightly larger final density at this momentum of about n(p = 2p ts ) ~ 0.5. We observe a 
good agreement of quantum and classical evolution for this momentum range as shown 
in Fig. |2|. The right figure shows the time dependent inverse slope parameter 

T(t,p) = -n(t,e p )[n(t,e p ) + l](dn/de)- 1 . (9) 

This parameter is constant for a Bose-Einstein distributed particle number and remains 
momentum dependent for classical thermal equilibrium. In Fig. |2] we plot the function 
T(t,p) for p low ~ and Phigh — 2p ts - Initially one observes a very different behavior of 
T(t,p) for the low and high momentum modes, indicating that the system is far from 
equilibrium. The quantum evolution approaches quantum thermal equilibrium with a 
momentum independent inverse slope T qm = 4.7 to high accuracy. In contrast, in the 
classical limit this slope parameter remains momentum dependent and the system relaxes 
towards classical thermal equilibrium, as exemplified in Fig. |l| 

I thank G. Aarts for collaboration on this work. 
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